RLE  COP?  AD-A152  616 


V  AFWAL~TR-84-30A5 


USE  OF  QUATERNIONS  IN  FLIGHT  MECHANICS 


WILBUR  L.  HANKEY 
L.  EARL  MILLER 
STEPHEN  J.  SCHERR 


Aerodynamics  and  Airframe  Branch  (AFWAL/FIMM) 
High  Speed  Aero  Performance  Branch  (AFWAL/FIMG) 
Aeromechanics  Division  (AFWAL/FIM) 


March  198A 


Final  Technical  Report 


Approved  for  public  release;  distribution  unlimited 


FLIGHT  DYNAMICS  LABORATORY 
AIR  FORCE  WRIGHT  AERONAUTICAL  LABORATORIES 
AIR  FORCE  SYSTEMS  COMMAND 
WRIGHT-PATTERSON  AIR  FORCE  BASE,  OHIO  A5A33 


85  '  03  28  055 


NOTICE 


Whan  Govarnmant  drawings,  specifications^  or  other  data  ara  used  for  any  purpose 
other  than  in  connection  with  a  definitely  related  Government  procurement  operation, 
the  United  States  Government  thereby  incurs  no  responsibility  nor  any  obligation 
whatsoevar/  and  the  fact  that  the  government  may  have  formulated,  furnished,  or  in 
any  way  supplied  the  said  drawings,  specifications,  or  other  data,  is  not  to  be  re~ 
garded  by  implication  or  otherwise  as  in  any  manner  licensing  the  holder  or  any 
other  person  or  corporation,  or  conveying  any  rights  or  permission  to  manufacture 
use,  or  sell  any  patented  invention  that  may  in  any  way  be  related  thereto. 

This  report  has  been  reviewed  by  the  Office  of  Public  Affairs  (ASD/PA)  and  is 
releasable  to  the  National  Technical  Information  Service  (NTIS) .  At  NTIS,  it  will 
be  available  to  the  general  public,  including  foreign  nations. 


This  technical  report  has  been  reviewed  and  is  approved  for  publication. 


JOSE^tTsTsHAiG  ^ 

Project  Engineer 


TOMMY  J.  MaJ,  USAF 

Chief,  Aer^ynamlcs  &  Airframe  Branch 
Aeromechanics  Division 


FOR  THE  COMMANDER 

RALPH  W.  HOLM,  Col,  USAF 
Chief,  Aeromechanics  Division 


* 


"Jf  your  address  has  changed,  if  you  wish  to  be  removed  from  our  mailing  list,  or 
if  the  addressee  is  no  longer  employed  by  your  organization  please  notify AFVJAL/FIMM  , 
w~PAFB,  OH  45433  to  help  us  maintain  a  current  mailing  list". 


Copies  of  this  report  should  not  be  returned  unless  return  is  required  by  security 
considerations ,  contractual  obligations,  or  notice  on  a  specific  document. 


UNCLASSIFIED _ _ 

•leUUlTY  CLA««l>IC*TION  0>  THU  »AOt 


TTTilToSTtlCUWrVCLAMIP  ICATION 

UNCUSSIFIED  _ 


at.  tieuMitV  CLAMirieATION  AUTHOAITV 


\0' 


REPORT  DOCUMENTATION  PAGE 

.  lib.  AiilAICTIVI  MAAKINOi 


2b.  OtClAMIbICATtON/DOWNOAADiNO  tCHIOUkI 


4.  MAFOAMINQ  OAOANIZATION  MAOMT  NUMMA(8> 


|3.  OI*TAI*UTION/AVAH.A»ltlTY  OA  AttAOAT 

Distribution  Unlimited 


t.  MONITOntNO  OnOANiZATION  AiAOAT  NUMIt  AISI 


AFWAL-TR-84-3045 


bt.  NAMI  OA  ABAAOAMINO  OAOANIZATION 

Flight  Dynamics  Laboratory 
AFWAL/FIMM  _ 


be.  ADDRtSt  (City,  SMI#  tub  Zlf  Cod#) 

Flight  Dynamics  Laboratory  (AFWAL/FIMM) 

AF  Wright  Aeronautical  Laboratories  (AFSC) 
Wright-Patterson  Air  Force  Base,  OH  A5A33 


b.  OAAICSSYMSOI.  |7*.NAMI  OA  MONITOAINO  OAOANIZATION 

afwttmM 


AFWAL/FIMM 


7b.  AOOASSS  rcily,  SmM  wid  ZIP  Cdd#l 


•t.  NAMI  OA  AUNDINQ/SAONSOAINO 
OAOANIZATION 


Be.  AOOAESS  ICily.  SMI*  and  ZIP  Codtl 


Sb.  OAAICI  SYMBOL  jS.  AAOCUAIMINT  INBTAUMENT  lOINTIAICATION  NUMIEA 

Of  tpplMdbl#)  I 


10.  BOUACi  OA  AUNDINO  NOS. 


AAOOAAM 
SLR  Ml  NT  NO. 


AAOJICT 

NO. 


WOAK  UNIT 
NO. 


2307 

1 1.  TITLE  (/ntJtd#  Stcurily  CMmiA##<io"I  UINUijAaairxcu  i 

USE  OF  QUATERNIONS  IN  FLIGHT  MECHANICS  » 


13.  AEASONAt  AUTHOAdI 

Wilbur  L.  Hankey,  L.  Earl  Miller.  Stephen  J.  Scherr  _ 

13a  TYAE  OA  AEAOAT  I3b.  TIME  COVEAEO  OATI  OA  AEAOAT  lYr.,  Mo.,  Day) 

Final  Technical  Report  t..i  go  ToSen  83  |  1984  March 


16.  SUAALBMENTAAV  NOTATION 


I  IB.  AAOE  COUNT 


C08ATI  CODES 


II.  SUBJECT  TIAMS  iContinu*  on  rwv*r§t  if  ntetuory  tnd  Idfntiry  by  bloc*  numbtri 


IS.  ASSTAACT  iConlinu*  on  nv«r««  if  ntetuary  and  idtntlfy  by  Woe*  number; 


Classical  flight  mechanics  employ  velocity  coordinate  or  body  axes  systems.  The 
latter  uses  Euler  angles  to  define  the  aircraft  orientation.  Singularities  and 
ambiguities  occur  whenever  the  flight  path  or  elevation  angle  Is  plus  or  minus  ninety 
degrees.  These  problems  are  circumvented  through  the  use  of  quaternions.  Two  problems 
are  addressed;  the  first  Is  a  nonsymmetric  body  and  the  second  is  a  symmetrical  configura¬ 
tion.  The  rotational  dynamics  and  the  solution  for  the  quaternions  and  Euler  angles  are 
determined  for  both  problems. 


ao.  OiSTAIBUTION/AVAILABItlTV  OP  ABSTRACT 
UNCLASBiPiBO/UNtIMITEO  IS  SAME  AS  APT.  C  OTIC  USB  AS  □ 


33a  name  op  ABSPONSISLB  INDIVIDUAL 

Wilbur  L.  Hankey 


31.  ABSTRACT  SICUAITY  CLAISIP ICATION 

UNCLASSIFIED 


23b  TILIAHONE  NUMSEA 
l/ncludt  Am#  Cod*; 

(513)  255-2455 


32e.  OPAICE  SYMBOL 

AFWAL/FIMM 


■  BITIQN  OP  1  JAN  n  IS  ObSOLlTI. 


msmsi 


FOREWORD 

This  report  mss  prepared  by  personnel  of  the  Flight  Mechanics 
Division  of  the  Flight  Dynamics  Laboratory,  Wright  Aeronautical  Labo¬ 
ratories  (AFWAL/FIM),  Wright-Patterson  Air  Force  Base,  Ohio.  This 
research  was  accomplished  under  Task  2307N603,  "Computational  Aerodynamics". 
The  effort  reported  here  covered  the  period  July  1982  to  September 
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1983.  Dr.  Wilbur  Hankey  and  Steve  Scherr  (FIMM)  accomplished  the 
original  work.  Or.  Earl  Miller  (FIM6)  starting  In  June  1983,  revised 
and  extended  the  results  of  Technical  Memorandum  AFWAL-TM-83-178,  "Use 
of  Quaternion  to  Numerically  Analyze  Aircraft  Spin".  This  TM  was  used 
as  the  basis  for  this  technical  report. 
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SECTION  I 
INTRODUCTION 


In  analyzing  aircraft  spin  and  other  flight  maneuvers  with  large 
rotation  rates  based  upon  first  principles*  the  non-linear  aerodynamic 
equations  (Navler-Stokes)  must  be  coupled  with  the  rigid  -body  dynamic 
equations  of  motion  .  In  both  systems  of  equations  It  Is  desirable  to 
use  a  body  axes  system  to  describe  the  forces  and  moments,  and  utilize  a 
surface-oriented  coordinate  system  for  obtaining  the  flowfleld  grid.  To 
solve  this  system  of  equations  It  Is  necessary  to  describe  the  position 
of  the  aircraft  with  respect  to  a  fixed  Inertial  coordinate  system.  The 

r, - 

classic  method  Is  to  use  the  Euler  angles  U-e-i|}  to  define  the  aircraft 
orientation.  However,  singularities  and  ambiguities  exist  when  the 
elevation  angle  Is  plus  or  minus  ninety  degrees.  This  difficulty  has  been 

overcome  In  the  field  of  Inertial  guidance  b^  the  use  of  quaternions  to 
describe  the  aircraft  position  In  space^^'^Tu  Is  the  purpose  of  this  re¬ 
port  to  Investigate  this  technique  for  numerically  solving  aircraft  spin 
problems.  Two  problems  are  addressed;  the  first  Is  a  nonsymmetrlc  body 
and  the  second  is^symmetrlc*-/^  The  rotational  dynamics  and  the  solution  for 


the  quaternions  and  Euler  angles  are  determined  for  both  problems. 
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SECTION  11 

V  DEFINITION  OF  QUATERNIONS 

) 

Hamilton  (circa  1840)  was  the  first  to  point  out  that  the  three 

■i 

Euler  angles  are  Inadequate  to  uniquely  define  the  orientation  of  a  body 
In  space  and  that  four  variables  are  required  to  resolve  the  predica¬ 
ment,  Hence  he  Invented  the  quaterni^ 


9^'A^quaterfilen-  Is  a  scalar  plus  a^ vectofTtotalling  four  elements. 


Q^t+v  :\ 

or  ,  ' 

Q  ■  (qj^,  Qj.  <ly  q^)  . 

Three  of  the  elements  describe  a  vector,  and  this  vector  defines  an  axis 
of  rotation.  The  fourth  element,  a  scalar,  defines  the  magnitude  of  a 
rotation  angle  about  the  vector  axis. 

The  definition  of  the  four  elements  of  the  quaternion  follows: 

-  cos  ^ 

q2  •  sin  ^  cos  a 

(2) 

q^  ■  sin  cos  6 
q^  -  sin  ^  cos  y 

where  a,B,Y  are  the  angles  between  the  rotation  axis  and  the  x,y,z  axis 
respectively.  The  amount  of  rotation  Is  the  angle  u.  The  rotations 
are  Illustrated  in  Figure  1* 

Observe  that  the  quaternions  are  subject  to  the  constraint  that 
Rj  +  qj  +  43  +  qj  •  1  (3) 

*  Figures  begin  on  page  23. 
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The  Interrelationship  with  the  Euler  trigies  Is  as  folipws: 


q.  •  C08  -I  cos  -f  cos  4  +  •!«  4  "f  2 


2 

1 


Qj  •  COS  -J  COB  -j  sin  2 


4  -  sin  #ln  ^  cos  ^ 


q^  ■  cos  -I  sin  |  cos  ^  +  min  ^  cos  ^  sin  ^ 


(4) 


sin  ■!  cos  ^  cos  "I  •  cos  sin  ^  sin 


The  relationships  for  the  Euler  angles  in  terms  of  quaternions  are 
sin  0  »  2(q^q3  -  qjq^) 
sln^l.  cos  0  «  2(q^q2  +  q3q^) 
cos  ♦  cos  0  ■  l-2q|  -  2q3 
sin  i|>  cos  0  ■  2(q^q^  4  qgq^) 


(5) 


cos  '!>  cos  0  ■  l-Zqj  -  2qJ 


The  four  quaternion  position  coordinates  are  computed  by  solving  four 
differential  rate  equations. 


e 

•*1 

0  -P  -q  -r 

‘’l 

^2 

1 

P  0  r  -q 

*>2 
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"  I 

^3 

q  -r  0  p 

“3 

• 

**4 

r  q  -p  0 

•^4 

(6) 


where  p,  q,  r  are  the  rotation  rates  In  roll,  pitch  and  yaw  respectively 
that  most  be  obtained  frtw  the  rigid  body  dynamic  equations.  These 
equations  are  well  behaved  since  p,q,  and  r  are  finite. 
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SECTION  III 

RIGID  BODY  DYNAMIC  EQUATIONS 

4 

The  rigid  body  dynamic  equations  used  to  compute  aircraft  spln^ 
with  six  degrees  of  freedom  contain  three  moment  equations  and  three 
force  equations.  The  selected  axes  system  for  the  differential  equation 
formulation  Is  a  body  coordinate  frame.  This  system  Is  selected  because 
the  aerodynamic  data  and  the  body  angular  rates  are  measured  relative  to 
It.  The  X  and  z  axes  form  a  plane  of  symmetry.  It  Is  convenient  to  align 
the  coordinate  system  with  the  principal  axes  since  all  products  of  Inerti 
are  zero  In  this  framework, 

Moment  Equations 

I„  P  +  (I  -I  )qr  •  l(u,VtW,p,q.ryp,S,b) 

X  *  y 

I  q  +  "  *l(u,v,w,p»q,r,p,S,c)  (7) 

y  X  z 

r  +  (ly-Ijj)pq  •  W(u,v,w,p,q,r,p,S,b) 

Force  Equations 

iD(u+qw-rv)-hng  sin  0«X(u,v,w,p,q,r,p  ,S) 

in(v+ru-pw)-iiig  sin  cos  e«Y(u,v,w,p,q,r,p,S)  (®) 

iB(^pv-qu)-mg  cos  cos  9"Z(u,v,w,p,q,r,p,S) 

Given  the  aerodynamic  characteristics,  along  with  the  mass  and 
Inertial  parameters,  the  above  six  dynamic  equations  are  solved  for 

p.q»r,u,v,w.  The  quaternions  are  obtained  by  solving  Equation  (6) 


as  a  subroutine  In  conjunction  with  the  simultaneous  solution  of  Equations 
(7)  and  (8). 


■«  *’r;  U  '».itf^.sr!»  -v^  3%  > 


The  hypothesis  at  this  point  It  that  the  numerical  solution  of  the 
quaternion  rate  Equation  (6),  Is  better  behaved  than  the  alternative 

procedure  for  resolving  the  Euler  angles.  The  differential  equations 
for  the  Euler  angles  are 


^  (q  aln  i  x  co*  a)  sec  G 
§  «  q  coa  i  •  T  tin  a 
a  ■  p  +  a  aln  © 


(9) 


Clearly,  difficulties  arise  when  e  reaches  plus  or  minus  90®.  The  rate 
for  the  a2<muth  angle  grows  without  bound,  thus  numerical  Integration  of 
Equation  (9)  Is  not  possible.  But  this  problem  1s  circumvented  by  first 
solving  for  the  quaternions  and  then  the  Euler  angles  from  Equations  (5). 

We  next  turn  our  attention  to  the  development  of  the  rotational 
solution  for  two  problems. 
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SECTION  IV 
MODEL  PROBLEM 


Nonaxi symmetric  Configuration 

The  first  problem  was  selected  to  demonstrate  the  essential  features 
of  manipulating  the  quaternions  and  to  obtain  the  required  numerical 
experience.  Observe  that  the  aerodynamic  moments,  Equations  (7),  are 
coupled  to  the  force  equations,  Equations  (8),  through  the  velocity  components 
For  a  situation  in  which  the  aerodynamic  torque  vanishes,  the  dynamic 
equations  are  uncoupled.  We  will  consider  such  a  case,  that  of  a  thin 
retangular  plate  which  is  illustrated  in  Figure  2(a).  In  this  case, 
the  moments  of  inertia  are  also  simplified  as  follows: 

I^  «  ly  +  I^  and  ly  >  (10) 

This  example  corresponds  to  a  configuration  where  all  products  of  inertia 
are  zero,  which  is  a  fair  approximation  for  most  aircraft.  For  the  case 
of  no  external  torque,  it  can  be  shown  that  rotation  about  the  y  vehicle 
axis  with  the  intermediate  value  for  the  moment  of  inertia  (ly)  is  unstable, 
while  rotation  about  the  other  two  axes  is  stable.  For  this  problem  the 
rotation  is  about  the  unstable  axis. 

Equation  (7)  now  reduces  to  the  following  simple  form: 


P  +  qr  ■  0 


q  -  pr  o  0 
•  2 

r  4  6  pq  ■  0 


(11) 


where 


e 


2 


ly  -  Ix 
ly  +  Ix 


<  1 


(12) 


Equation  (11)  combined  with  Equation  (6)  will  describe  the  angular  motion 
of  the  rigid  body  given  a  set  of  initial  conditions.  (These  values  can 


then  be  Inserted  Into  Equation  (8)  to  find  the  translational  components). 


Initial  conditions  for  the  unstable  situation  were  selected. 


p(0)  -  cO 
q(0)  •  n 
r(0)  «  0  . 

Q(0)  -  (1.0, 0.0) 


Axl symmetric  Configuration 

The  shape  for  this  problem  is  a  thin  circular  disk  with  the  x  and  y 
axes  In  the  circular  section  and  i  perpendicular  to  these  axes.  The 
configuration  Is  Illustrated  In  Figure  2(b).  All  products  of  Inertia 
are  zero.  The  moments  of  Inertia  are 


•  ly  •  mR^ 
7  mR^ 


(14) 
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SECTION  V 
NUMERICAL  RESULTS 


Nonaxlsyrmietrlc  Configuration 

The  system  of  equations  (Equations  (6)  and  (11))  were  solved  for 
parametric  values  of  fl  ■  10  rad/sec.  t  •  0.01  and  d  ■  0.7. 

A  Runge-Kutta  fixed  time  step  Integration  technique  was  used. 

Figures  3»  4i  5  present  the  roll*  pitch  and  yaw  rates  for  several 
revolutions  as  computed  from  Equation  (11).  Figures  6  through  9  present 
the  history  of  the  quaternion  position  coordinates  as  determined  from 
Equation  (6)  using  the  values  of  p.  q,  and  r  determined  from  the  solution 
of  Equation  (11).  The  time  step  size  used  for  the  Integration  was  0.001  sec 

which  was  selected  after  numerical  experimentation  produced  repeatable 
results.  The  Euler  angles  were  then  computed  from  Equations  (5). 

Figures  10  through  12  present  the  history  of  the  Euler  angles.  No 

difficulty  was  encountered  In  numerically  Integrating  Equations  (6)  and 
(11).  The  reason  for  this  Is  that  the  quaternions  are  obtained  from  Equations 

(6)  which  are  well  behaved.  Also,  they  are  bounded  by  plus  and  minus 
one.  Thus  Equations  (11)  are  well  behaved  and  no  difficulty  Is  obtained 
In  Integrating  them.  The  difficulty  Is  In  attempting  to  Integrate 
Equations  (9).  But  this  Is  unnecessary  since  the  Euler  angles  can  be 

obtained  directly  from  Equations  (5). 

The  solutions  for  the  Euler  angles  i|>  and  (t*  undergo  rapid 
transitions.  This  Is  a  result  of  two  conditions  that  take  place.  From 
Equation  (5)  ,  6  approaches  plus  or  minus  90”  and  the  combination  of  the 
quaternions  on  the  right  side  changes  sign.  Thus  the  right  side  changes 
rapidly  from  large  positive  values  to  large  negative  values  and  vice 
versa. 
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A  further  observation  Is  that  the  rotational  rates  are  periodic 
functions  of  time.  But  the  Euler  angles  are  not. 

An  analytical  solution  to  this  problem  exists  (see  Appendix  A). 
Excellent  agreement  was  obtained  between  the  numerical  and  analytical 
values  thereby  confirming  the  numerical  approach. 


Axlsymmetrlc  Configuration 

The  system  of  equations*  Equations  (6)  and  (16)  were  solved  for  a 
parametric  yaw  rate  of  n  *  10  rad/sec.  It  Is  shown  In  Appendix  B  and 
can  be  verified  numerically  that  the  solution  for  the  body  angular  rates 


p  ■  0 

(19) 

r  “  fi 

Figures  13  through  ]6  ^^f"®  histories  of  the  quaternions.  Figures 

17  through  19  are  the  Euler  angles. 

Two  significant  results  are  established  In  Appendix  B.  The  first 
Is  that  the  configuration  remains  In  the  same  vertical  plane  that  it 

started  In  whenever  the  bank  angle  Is  zero.  The  second  result  Is  that 
the  attitude  does  not  change  throughout  the  trajectory.  These  results 

agree  with  the  observed  trajectory  of  a  Frisbee  when  the  Initial  bank 
angle  is  zero.  The  numerical  solution  agrees  with  the  analytical  solution 
of  Appendix  B. 
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SECTION  VI 
CONCLUSIONS 

Experience  has  shown  that  numerical  difficulties  arise  In  using  the 
Euler  angles  to  describe  the  orientation  of  an  aircraft  during  large 
changes  in  the  Euler  elevation  angle.  The  angular  rates  for  the  heading  and 
bank  angles  are  undefined  whenever  the  elevation  angle  reaches  plus  or  minus 
ninety  degrees.  The  numerical  difficulties  may  be  completely  overcome 
through  the  use  of  quaternions  In  place  of  the  Euler  angles.  The  Euler 
angles  are  then  easily  computed  from  the  quaternion  solution. 

Two  problems  were  solved  numerically  and  analytically.  The  results 
for  the  two  approaches  did  agree  and  there  was  no  difficulty  In  obtaining 
numerical  results.  The  first  problem  was  an  example  of  a  nonaxi symmetric 
configuration.  The  rotation  was  about  the  unstable  principal  axis.  The 
body  rotational  rates  were  periodic,  however,  the  quaternions  and  Euler 
angles  were  not. 

The  second  problem  represented  an  axisymmetric  configuration.  Again 
no  numerical  difficulties  were  obtained  when  the  quaternions  were  used. 

The  rotation  was  about  the  yaw  axis.  Two  significant  results  obtained 
were  that  if  the  initial  bank  angle  was  zero,  then  the  trajectory  remained 
in  the  initial  vertical  plane.  The  second  result  was  that  the  attitude  did 
not  change  throughout  the  trajectory.  The  use  of  quaternions  to  solve  flight 
mechanics  problems  Is  therefore  highly  recommended  for  future  Investigations. 
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APPENDIX  A 


ANALYTIC  SOLUTIONt  NONAXISYMMETRIC  CASE 
An  analytic  aolutlon  axlata  for  tha  nodal  problan  aalaetad. 


p  +  qr  -  0 

p(0)  •  cn 

q  -  pr  ■  0 

q(0)  -  fl 

(A-1) 

k..  „ 

T  +  fi^pq  ■  0 

r(0)  -  0 

Regroup  these  three  moment  equations  as  follows: 

p  ■  -qr  ■  -  qq/p  or 

^  (p^+q^)  ■  0 

L 

q  ■  pr  »  -  rr/6  q  or 

(rVfi^q^)  -  0 

(A-2) 

r  ■  -6  pq  ■  6  p?/r  or 

IF  («V-t2)  -  0 

These  three  equations  may  be  integrated  readily  and  evaluated  using  the 

initial  conditions.  Hence, 
+  <1^  • 

+  rW  . 

-  c'tl^ 


(A-3) 


These  expressions  are  energy  constraint  relationships,  however,  only  two 
independent  equations  are  evident  with  the  third  one  being  redundant.  There¬ 
fore,  one  additional  differential  equation  nust  be  solved  and  the  two  energy 
constraints  used  to  obtain  the  complete  solution. 

Differentiate  Equation  (11)  with  respect  to  time. 


"  •  •  2  2  2 
p  ■  -qr  -  rq  -  -p  (r  -  6  q  ) 

.  2  2  2 
q  »  rp  +  pr  •  -q  (r  +  6  p  ) 

,2*  .2*  .2  ,22. 

r  ■  -  6  pq  -  6  qp  ■  6  r  (q  -p  ) 


(A-4) 


The  energy  constraints  of  Equation  (A-3)  may  be  inserted  into  one  of  the  above 
equations  to  obtain  a  single  differential  equation  with  one  dependent  variable, 
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Wc  vill  select  the  nlddls  equation. 

q  +  ad^Id+c^/j)  ft*  -  q^lq  -  0  (A-5) 
This  equation  is  called  "Ouffinge  Equation"  end  possesace  an  analytic 
solution.  The  solution  is  obtained  as  follows: 


Observe  that 


q 


Hence  Equation  (A-6)  becomes: 


d  q*  -  -  46*  l(l+c*/2)  ft*  -  q*lq  dq 


Equation  (A-7)  may  be  Immediately  integrated  once 
q*  -  6*q*  Iq*  -  (2  +c*)ft*l  + 


(A-6) 


(A-7) 


(A-8) 


The  initial  conditions  may  be  used  to  evaluate 


q(0)  -  p(0)  r(0)  -  0 
q(0)  ■  ft 

Hence 

-  6*(l+e*)ft^ 

Equation  (A-8)  then  becomes 

q*  -  6*(n*  -  q*)  (f.*(l+c*)-q*J 


(A-9) 


(A-10) 


(A-11) 


’  Solving  for  q  gives 

q  -  +  6/2*-q*  V^*(l+c*)-q* 

The  sign  Is  determined  as  follows:  Initially 
q(0)«n  >  0 
q(0)  -  0 

q(0)  -  -  q(0)  6*P*(0)  <  0 


(A-12) 


14 


(A-13) 


Thus  In  the  neighborhood  of  t«0,  q  Is  a  convex  function.  It  follows 
that  for  t>0,  q<0  and  therefore  the  negative  sign  Is  selected. 

A  convenient  transformation  Is  next  employed.  Let 

Bln  X  -  ^  (A-l 

Differentiating  gives 


1  ^  -  fl. 

cos  X  ^  -  fl 


Since 


X(0)  -  e/2 


and  q  decreases  with  time 
X(t)  <  Tr/2 

Equation  (A-16)  can  be  written  as 


(A-17 


(A-ia 


q  .  ‘  X 


(A-19 


Equating  to  Equation  (A-14)  gives 


i  .  .  M  JW  Bln"x 
k 


where 


k"^  -  1+e^ 


Equation  (A-20)  can  be  Integrated 


Rewrite  the  right  side 


A  table  for  the  solution  follows  for  a  k  «  0.99995 


X 

%/a~ 

S£L-t 

k 

90® 

5,991 

1.000 

0 

85® 

3.126 

.996 

2.865 

o 

O 

00 

2.436 

.985 

3.555 

o 

e 

1.735 

.940 

4.256 

60® 

1.317 

.866 

4.674 

50® 

1.011 

.766 

4.980 

40® 

.763 

.643 

5.228 

o 

o 

.549 

.500 

5.442 

o 

O 

CM 

.356 

.342 

5.635 

10® 

.175 

.174 

5.816 

0® 

0 

0 

5.991 

The  relations  for  p  and  r  are  determined  in  a  manner  similar  to  the 
way  that  q  was  determined.  Initially 

P(0)  -  0 

p(0)  •  p(0)6^q^(0)  ^  ®  (A-27) 

^(0)  <  0 

In  the  neighborhood  of  t*0,  p  Is  a  concave  function  and  r  Is  a  decreasing 
function.  From  Equation  (A-3) 


Since  p  Increases 
.  fi  jm 


P  ■  ^  vl-it^ain^A 


(A-29) 


From  Equation  A-3 

ira-  ±  S 

n  ±  <:-•>»  A 

From  the  inU^el  bourwiary  condition  and  AoVcCt  it  follows  that 


(A-30) 


r  ■  -  6fl  cos  X 


(A-31) 


The  following  table  presents  the  numerical  values  obtained  by 
numerical  Integration. 


Cic-c.) 

? 

C'T^.ti/s^C  } 

^TrAeZ/seC  J 

■  '  •  ~  . .  " 

. 

0 

0.10 

10.00 

0 

.409 

0.88 

9.96 

-0.61 

.508 

1.74 

9.85 

-1.22 

.608 

3.42 

9.40 

-2.39 

.668 

5.01 

8.66 

-3.50 

.711 

6.41 

7.68 

-4.49 

.747 

7.66 

6.42 

-5.36 

.777 

8.65 

5.02 

-6.05 

.805 

9.40 

3.42 

-6.58 

.831 

9.85 

1.73 

-6.90 

.856 

10.00 

0 

-7.00 

*  •  . 

•  f- 


1  \  • 


':'*r''.V7*.Tn.\x.  •  n.-^  m-,  k  ■  .cTg.T  iCF»  ■ 


The  agreement  between  the  theoretical  values  for  the  body  rotational 

rates  and  those  obtained  by  numerical  Integration  of  the  differential 
equations  is  very  good.  This  can  be  observed  by  comparing  the  tabular 
results  with  Figures  3,4,  and  5,  Observe  that  the  “rate  period”  Is 


^Ic 

Period  -  F(iT/2,k)  -  3.423  aec 


(A-32) 


and 


Note  that  the  range  may  be  extended  as  follows  for  any  Integer  j; 
PCjtt+X.k)  -  2jF(7r/2,k)  +  F(X,k)  (A-33) 

F(-X,k)  -  -  F(X,k)  (A-34) 


t 


K-- 

r 
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C 


t 


f' 


J*'. 

s\ 


V 


t 
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APPENDIX  B 

ANALYTIC  SOLUTION:  AXISYMMETRIC  CASE 


The  body  angular  rates  satisfy 

p  4  qr  •  0  P(^)  "  0 

q  -  rp  -  0  <l(0)  ■  0 

f  ■  0  ^(0)  • 

From  the  latter  equation 
r  ■  fl 

Combining  the  first  two  equations  gives 

•  2 

P  “  -qr  -  -r  P 

The  boundary  conditions  are 

p(0)  -  p(0)  «  0 

The  solution  Is 

p«0 

In  a  similar  way,  It  can  be  concluded  that 
q*0 

Substitution  of  p,q,  and  Into  Equation  5  gives 

’l  ■  -  I  ■^’•4 
’2  ■  I  '•I3 
’3  ■  ■  5  '■‘I2 
’4  ■  I  '■’l 


(B-1) 

(B-2) 

(B-3) 

(B-4) 

(B-5) 

(B-6) 

(e-7) 


Combining  the  first  and  last  equations  gives 
**  /I 

**1  "  (B-8) 

Analytic  Integration  using  the  boundary  conditions  gives 

%  tit  ,  , 

qi  •  cos  j-  cos  Y~  (g.g) 

can  be  obtained  directly  from  the  first  differential  equation 


2  « 


cos 


®0 

T- 


sin 


(It 

7” 


(B-10) 


-g.* 
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Combining  the  second  and  third  differential  equations  results  in  a 
second-'Order  differential  equation.  Integrating  and  satisfying  the 
boundary  conditions  gives 


q,  ■  sin  T—  sin 


0 

■  sin  cos 


The  Euler  angles  are  obtained  from  Equation  (14) 
sin  0  -  2(qj^qj  -  ®0 

Bln  ♦  -  2(qj^q2  **■  ®o  ® 

2  2 

cos  ♦  ■  (l'2q2  -  2q2)  sec  ©  •  cos  0^  sec  0 

sin  tp  ■  2(q2q2  +  **1**4^  ®  "  ® 

2  2 

cos  ip  ■  (l“2q2  -  2q^)  sec  0  ■  cos  0q  cob  fit  sec  0 

Let  the  inertial  coordinate  axes  be  defined  as  follows:  The  : 
axes  are  In  the  vertical  plane  of  the  Initial  motion.  The  y^ 
perpendicular  to  these  axes.  The  kinematic  relation  for  y^  1i 
the  quaternions  and  velocity  components  Is 

yg  »  2u(qj^q^  +  q2q3)  +  v(l-2q2  -  2qJ)  +  2w(q2q^  -  q^q^) 
■  u  sin  fit  +  V  cos  fit 

Let  Vp  denote  the  component  of  the  velocity  vector  measured  1i 
plane.  Then  the  components  u  and  v  satisfy 

u  ■  Vp  cos  fit 


(B-11) 

(B-12) 

^  and  Z£ 
axis  Is 
I  terms  of 

(B-13) 

I  the  x-y 


V  — sin  fit 


(B-14) 


These  two  equations  follow  from  the  result  that  the  x  and  y  axes  are 
rotating  at  a  constant  rate  Q.  Substitution  of  u  and  v  Intoy^  gives. 
y^^O. .  Therefore  y^  is  t:::ual  to  Its  Initial  result  and  does  not  move  out 
of  the  plane. 

Since  p  and  q  are  zero,  the  angular  velocity  of  the  x,y,z 
coordinate  frame  Is  u  where 


w  *  re_ 

The  rate  of  change  of  1s  determined  from 
4-e 

dt  ^  -e  -a 


(B-15) 


(B-16) 


Therefore  which  lies  along  the  z-axis  does  not  change. 
Consequently  the  orientation  of  the  z-axIs  does  not  change.  The 
conclusion  Is  that  the  attitude  does  not  change. 


Figure  8 
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